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Abstract 

We calculate the integration ideal of annihilating differential operators 
of the non-normalized Fisher-Bingham distribution and show that the 
ideal agrees with the set of operators for the Fisher-Bingham integral 
given in [9]. They conjectured that the set generates a holonomic ideal 
and we prove their conjecture. 

1 Introduction 

The Fisher-Bingham distribution is a probability distribution on the n-dimensional 
sphere S n (r) with the radius r defined by 

' exp(t T xt + yt)\dt\. (1.1) 



F(x,y,r) 



Here, the variable x is an (n + 1) X (n + 1) symmetric matrix whose 
component is Xij when i = j and Xij/2 when i ^ j. The variable y (resp. t) is 
a row (resp. column) vector of length n + 1, and the measure \dt\ is the Haar 
measure on S n (r). The function F(x,y,r) is the normalizing constant defined 

by 

F(x,y,r) = J exp I ^ XijUtj ■ + ^ i^U J \dt\. (1.2) 



l<i<j<n+l 



The integral (|1.2|) is referred to as the Fisher-Bingham integral on the sphere 
S n (r). 

The Fisher-Bingham distribution is used in directional statistics. Kent stud- 
ied estimations, hypothesis testings and confidence regions with respect to the 
Fisher-Bingham distribution on the 2-dimensional sphere [3], and in the book 
by Mardia and Jupp on directional statistics [H chapter 9] , a definition of the 
Fisher-Bingham distribution having the same form as (|1.1|) and a relation with 
an another probability distribution on the sphere are explained. 
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We are interested in estimating the value of parameters x and y which max- 
imizes the likelihood function 



f( x > v) = w II ex p(^ xt i + y*i) 

for given data t±, ■ ■ ■ ,tpf 6 S n . This problem is equivalent to estimating the 
value of parameters x and y which minimizes the function 

F(x, y,r)exp - ^ S^x.^ - S i y l 

\ l<i<j<n J 

for given {Sy}i<,'<j< n , { -S'-i } 1 < i < n C R. There are several approaches to solving 
this problem. Among them, the holonomic gradient descent proposed in [9 
enables us to estimate the value by utilizing linear partial differential operators 
with polynomial coefficients which annihilate the Fisher-Bingham integral (ll.2[) 
and generate a holonomic ideal. Let Dd be the ring of differential operators 
Dd = C{z\, . . . , Zd, d\, . . . , dd)- A left ideal in Dd is called a holonomic ideal 
when the characteristic ideal m^ ^(I) generated by the principal symbols of 
/ in C[zi, . . . , Zd, £1, • • • ,£d] has the Krull dimension d. See, e.g., p 31, 
Definition 1.4.8], [BJ, and their references for details. 

In [5], it is shown that the Fisher-Bingham integral F(x, y, r) is annihilated 
by the following linear partial differential operators. 

d Xij -d Vi d Vj {i<j), 

n+l 

i=l 

+ E ( x kjd Xik ~ x ik d Xjk ) + yjd Vz - y l d V] (i < j,x k t = xtk), 
rd r - 2 ^ ■•'■j'lr - E V id Vi ~ n - 

i<j i 

They also show that (|1.3[) generates a holonomic ideal in the cases n = 1 and 
n = 2 by a calculation on a computer, and conjecture that it holds for any n. 
We will prove this conjecture. 

In order to state the main result of this paper precisely, let us explain the 
notion of the integration ideal. For a holonomic ideal I in Dd, the left ideal 
(/ + dd'+iDd + ■ ■ ■ + ddDd) H Dd' in Dd' is called the integration ideal and it 
is known that the integration ideal is a holonomic ideal in D^ (see, e.g., [3J 
Chapter 1], §5.5]). 

In the present paper, we show that (|1.3[) generates the integration ideal of 
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the annihilating ideal 



Ann exp ^ 

\ \l<i<i<n+l » =1 / / 

Here, {z±, . . . , Zd>} = {xij,yk\l < i < J < n + 1, 1 < fc < n + 1} and 
{^d'+ij • • • j Zd} — {ti, ■ • ■ , t-n+i}- As its corollary, we show that (ll.3[) gener- 
ates a holonomic ideal for any n and prove the conjecture in 9' . Oaku gave 
an algorithm for computing the integration ideal in [TJ. The proof for n = 1, 2 
are done by applying this algorithm on a computer. We apply this algorithm 
for a general natural number n, for which the steps of the algorithm cannot 
necessarily be applied, and so some propositions are necessary. 

In section [21 we consider the holonomic ideal annihilating the Haar measure 
on S n (r). In section [31 we give generators of the holonomic ideal which anni- 
hilates the integrand of the Fisher-Bingham integral. In section [H we compute 
the integration ideal of the holonomic ideal which is given in section [3l and 
prove the main theorem of this paper. 



2 The Haar measure on S n (r) 

The Riemannian metric on the n-dimensional sphere with radius r is constructed 
by the pullback of the standard metric on the (n + l)-dimcnsional Euclidean 
space R™ +1 along the embedding map. The metric defines a probability measure 
on S n (r), which is called the Haar measure and denoted by \dt\. We define a 
distribution \i r with a parameter r > as 

(Mr, </>(*)) := / 

Here, (pit) is a test function. 

Let D = C(x,y,r,t,d x ,dy,d r ,dt) be the ring of differential operators with 
polynomial coefficients. For a given distribution F, we denote by Ann(F) the 
set of the operators in D which annihilate F. 

Lemma 1. A left ideal I in D generated by following differential operators is 
a subset of Ann(/Lt r ). 

d Xij (1 < i < j < n + 1), d v .(l<i<n + l), tf + 

n+l 

Udtj - tjd u (1 < i < j < n + 1), rd r + ^ Ud, + 1 

i=l 

For computing the integration ideal, the following proposition is important. 

Proposition 1. The left ideal I in D is a holonomic ideal. 

This proposition may be well known, however we could not find a proof in 
the literature. Therefore, we present a proof here. 



i L n+1 ' i 

(2.1) 
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Proof. By the fundamental theorem of algebraic analysis (see, e.g., [5J p30.Theoreml.4.5]), 
it is enough to show that the dimension of the characteristic ideal in( 0)e )(J) is 
not more than the number of variables N := n(n + l)/2 + In + 1. 

We can find the operators r 2 dt h + tkrd r — tk (1 < k < n + 1) in / as follows. 

t n+1 (t n+1 d tn+1 +■■■ + hd tl + rd r + 1) - d tn+1 (t 2 n+1 + ---+t\-r 2 ) 

n 

= -^U(Ud tn+1 - t n+ id ti ) + r 2 d tn+1 + t n+1 rd r - t n+1 , 

i=l 

tk{t n +idt„ +1 H 1- tid tl + rd r + 1) - t n+ i(t k d tn+l - t n+ id tk ) 

k—1 n 
i=l i=fc+l 

+r 2 9 tfe + i fe r<9 r - t fc (1 < k < n) 
Then, the characteristic ideal in( 0je )(J) contains the polynomials 

£*„(!<<< j<n + l), S Vi (l<i<n + l), t 2 n+1 + ■ ■ ■ + t\ - r 2 , 
Htt ~ (1 < i < 3 < n + 1), r 2 i u + Ur£ r {l < i < n+ 1). 

Let /' be the ideal in the polynomial ring C[x,y,r,t,£ x ,£ y ,£ r ,£t] generated by 
these polynomials. Then, we have I' C in(o, e )(-0- Since dim/' > dimin(o e )(J), 
it is enough to show that dim I' < N. 

Consider the graded reverse lexicographic order satisfying 

&„+i > — y&ytxyZyy&y t n +i > — >- h> x>- r. 

Since the degree of the Hilbert polynomial of an ideal in the polynomial ring 
equals that of the initial ideal with respect to the graded order of the ideal 
(see, e.g., [1] p448, Proposition 4]), the dimension of I' is equal to that of 
the initial ideal LT^(J') with respect to this order. The initial ideal LT^(Z') 
contains the monomials , £ ai , ii£t., , , £ 2 +1 . Let I" be the ideal generated 
by these monomials. Analogously, we can show that it suffices to prove that the 
dimension of /" is not more than N . 

Computing the irreducible decomposition of the algebraic variety defined by 
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I" as 

V(UnZy k M ti ,r 2 Z tk ,t 2 n+1 ;l <k<l<n + l,l <i < j <n+l) 

n+l 

= ^,^,tn+i;i<Kj<n + i)n f) V(ut tj )n f) V(i*Z u ) 

\<i<j<n+l i=l 
n+l 

= %,4..Wi;i< % < j <n + i)n (J V(*i, . . . ,fi-i,£t i+1 , • • • )6 n +i) 

i=i 

n(F(r)UF(e tl) ...,( (n+1 )) 

/n+l \ 



vfe=l 



n(^(r)U7(f tll ..,U) 



'n+l 



= f IJ V r ($ Xfcl ,$ w ,r,t T ,+i ) *i,...,ti_ 1 ,^t <+ i,.-.,6„+i;l < fc < ' < n + !)J 

U V(e Xw ,^ fc ,t n +i,ti,...,t i _i ) Ct 1 ,...,6„ +1 ;l< fc < J < n + l)^ , 
we conclude that the dimension of I" is exactly N. □ 

3 Holonomic ideal annihilating exp(p)/x r 

Let g(x,y,t) be the polynomial Z)i<,-<j< n+1 scytitj + 2/»*«- We can S et a 

holonomic ideal annihilating the distribution exp(g(x,y,t))fj, r by the following 
lemma. 

Lemma 2. Consider the ring of differential operators with polynomial coeffi- 
cients C(xi, . . . , X n , di, . . . , d n ). Let u be a distribution and suppose that I C 
Ann(u) is a holonomic ideal. Let f be a polynomial and /, := df jdxi. Then, 
the left ideal J generated by 

;d Xl ,...,d x J e 1} 

is a holonomic ideal such that J C Ann(e-'w) 

For a proof of this lemma, we refer to [5] . It follows from this lemma that the 
left ideal J in D generated by the following differential operators is a holonomic 
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ideal and included in Ann(exp(g)/v). 

d Xij - tit j (l<i<j<n + l), 
d Vi -U (l<i<n+l), 

n+l n+l 
- J! ^' fc * fc ~~ ~~ ~~ " Xiktk ~ Xliti ~ 

k=1 s fc=1 (3-1) 

(l<«<j<n+l), 

*1 + ' ' ' + *n+l — r "i 

n+l / n+l \ 

rd r + 1 + *i I ~ ^ ^t*** ~ ^M** ~ f» J 

i=l V fc=l / 

In fact, we will show that the ideal J is generated by the differential operators 



U - d Vi (1 < i < n + 1), d Xij - d Vi d Vj (1 < i < j < n + 1), 

n+l 

+ in,, <L . - x ik d xjk ) + y 3 d y% - yid Vj + d u d Vj - d tj d m ( 3 - 2 ) 
(1 < i < j < n + 1, a;fci = Zft), 

n+l n+l 

rd r - 2 .r,.<7.. - ^ yid Vi - n + ^ d u d Vi 

i<j i—1 i—1 

To prove this statement, we prepare the following lemma. 
Lemma 3. 

t a = d° mod D{ti-d m ;l<i<n + l} (3.3) 

Proof. When a = e», the equation (|3.3[) obviously holds. Let us assume that 
(|3.3p holds for the case of a — ej. Then, we have 

= Ud^-^ mod D{ti-d n ;l<i<n + l} 

= d^~^U = df-^fc -dy i ) + df-^dyt 

= d^- ei) d y% mod D{U - d Vl ; 1 < i < n + 1} 

Hence, (|3.3p holds for a. Therefore, the equation Q3.3[) holds for any a. □ 

Finally, we prove the following lemma. 
Lemma 4. The differential operators p.2[) generates J. 
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Proof. Let K be the left ideal generated by p.2[) . First, let us show J C K. 
The equation 

9x w - Utj = d Xij - d Vi d Vj mod D{U - d Vi ; 1 < i < n + 1} (3.4) 

gives the inclusion — fjtj e K. 

The inclusion — t{ G if is obvious. The inclusion ti(dt j — x jktk — 

x jj l j ~Vi)~ tj{du ~ Sfeii ^fe*ft - ^lii* - Hi) € AT follows from 

n+1 n+1 



*i(^tj _ X! X J fc ** ~~ ~~ y ^ ~ tjiPu - ^ Xiktk ~ Xi%ti ~ 

k=l k=l 
n+1 n+1 

^ x ik^k^j ^ x jkXk^i ^jj^i^j x ii^i^j Vj^i ~t~ Ui^j ti$tj ^jdt i 

k=l k=l 

n+1 

^ yC^fe^j x jk^i)^k x jj^i^j yi^j ^i^tj ^jdf i 



k=l 
n+1 



— ^ ] { x ikdyj x jkd Vi )dy k + x jj)^Vi^Vj 



k=l 



+yid yj - + d^c^ - <9 w <9 ti mod D-j^ - 9^; 1 < i < n. + 1} 



n+1 



— ^^(^ife&jfc x jkd Xik ) + (xu x jj)d Xi 



fc=l 



+2/i<9 % - y 3 dy z + d Vi d tj - dy 3 d ti mod D{d Xij - d yi d y] ;l < i < j < n + 



x 



Since 



+Uid Vj - y 3 d Vt + d Vi d t] - d Vj d ti 



4 + ---+tl +1 -r 2 



dy 1 + --- + dy n+1 2 - r 2 mod D{U - d Vi ;l < i < n + 1} 



n+1 



^2 dx " ~ r2 mod D {&*v ~ 9y z d Vj ;l < i < j < n + 1}, 



i=l 



we have Y^=i ®xu ~ r2 e ^ v • 

The inclusion rd r + 1 + 22™=/ (ft, - ^Zfc=i x iktk - XuU - yij € K follows 
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from 

n+1 / n+1 

rd r + 1 + ^2 *» I du - £ifc*fc - aJiiti - 

i=l V fc=l 



n+1 

ra r + 1 + ^ 


n+1 n+1 


n+1 


n+1 


i=l 


i=l k=l 


i=l 


1=1 


n+1 


n+1 n+1 


n+1 


n+1 


r<9 r - n + d ti tj 

i=l 


^ ^ ^ ^ %ikt{tk 
i=l k=l 


1=1 


- 51 

1=1 



n+1 n+1 n+1 n+1 n+1 

= rd r ~ n + ^ ®U 9 vi - 2J X! Xtkdx ^ ~ 5Z ^"^m ~ 

2—1 i— 1 fc=l i= 1 i= 1 

mod D{U - d Vi , d Xij - d yi d Vj ; 1 < % < j < n + 1} 

n+1 n+1 

= rd r - n + ^ d tl d Vl - 2} r, ,il r - 

Z=l 2<j Z=l 

Therefore, we have J C K. 

Next, let us show the opposite inclusion K C J. The inclusion - 9 9i 6 J 
is obvious. The inclusion — d Vi d yj G J follows from equation p.4|) . Other 
generators of if are also in J because of the above equivalence relation. □ 

4 The Fisher— Bingham Integral 

Let D' be the ring of differential operators with polynomial coefficients C{x, y, r, d x ,d y , d r ). 
The left ideal J' := D' n ( J + {9 tl , . . . , 9t n+1 } ■ D) in £)' is the integration ideal 
of J. The Fisher-Bingham integral (|1.2p can be written as 

F(x, y, r) = {e^'V^Hr, 1} = f exp(g(x, y, t))fi r dt. 

JR"+1 

Hence, the operators in J' annihilate F(x, y, r). It is known that the integration 
ideal of a holonomic ideal is also a holonomic ideal (see, e.g., [3J 2, chapterl]). 
Therefore, if we obtain a set of generators of J', then this set generates a holo- 
nomic ideal. In this section, we compute a set of generators of J'. As the first 
step, we prove the following lemma. 

Lemma 5. Let P be an arbitrary differential operator in (|3.2p ," then we have 

t a P = dy'P mod D{U - d Vi ; 1 < i < n + 1}. 
Proof. For simplicity, put Qij = x ij d Xu +2(x jj -x li )d Xlj -x lj d Xjj +J2k^ l: j ( x kj9x ik ~ x lk d Xjk ) 
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and R = rd r — 2 Y]j<rj Xijd Xij — n. The following equations prove the lemma. 

'" (Qij + Vjdy - !!.<>„ + d ti d Uj - ih d y ) 
= (Qij + Vjdy, ~ Vidyj + 9 u d Vj - d tj d yt ) t a ~ ^dy.d^^ + aj%i9 t (Q ~ ej) 
= (Qij + Vjdy, ~ Vi d Vi + <>' 9 v " dt^) <)\] 

-ctidy^ 01 -**) + ajdy^-^ mod D{U - d Vi ; 1 < i < n + 1} 
= d y (Qij + Vj d y, ~ vAi + d u 9 yj ~ d t} d Vi ) , 



n+1 n+1 



h^-E^+E^a 



n+1 n-\-l \ n+1 



= [r - E Vid yi + E 9 *A Y a ~ E «*M a ~ ei) 

V i=l 1=1 / i=l 

/ n+1 n+1 \ n+1 

= \ R ~ E + E 9 *A ) d v ~ E ««%4 Q_e,) mod - ; l < i < n + 1} 



\ i=l i=l / 

(n+1 n+1 
i=l i=l 



□ 



Theorem 1. TTie integration ideal J' is generated by the differential operators 
m (fj~3l). 



Proof. Let F and F' be the set consisting of the differential operators Q3.2p and 
(|1.3[) respectively. The inclusion D' ■ F' C J' is obvious. We need to show the 
opposite inclusion D' ■ F' D J'. If a differential operator P is contained in J', 
then P can be written as 

P = E QvPi + E ^ ^ ( p < e F, Qi e D, Rj e D), 

i j 

from the definition of J'. Without loss of generality, we can assume that no 
term of Qi contains dt- Note that 

t a Pi ee d™ Pi mod D{t k - d Vk ; 1 < k < n + 1}; 

then, P can be written as 

p = E Q'^+E 9 *^+E ^C**-^) e F,Q' t e D', R j eD,S k e D). 

i j k 

Since all differential operators in F except U —d Vi have the form P'+X), ^t* U[ (P' G 
F 1 ,U[ e D'), P can be written as 

P = E Q^'+E 5 *^+E ( p e F > Qi e D 'i P J e D > S * e 
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Moving some terms to the left-hand side, we obtain 

Q'i p i-Yl S k{tk-d yk ) = ]T d t,Rj (Pi G F', Q\ G D', R 3 eD,S k e D) 

i k j 

Without loss of generality, if we assume that no term of Sk contains dt, then 
the left-hand side of the equation does not contain d t - Expanding both sides 
and comparing the coefficients, we get ^ . d tj Rj = 0, in other words, we obtain 

i k 

The right-hand side of this equation is included in the left ideal D ■ {U — d Vi |1 < 
i < n + 1} in D. Let the weight of ti be 1 and that of other variables be 
0, and consider a term order -< with this weight. The Grobner basis of D ■ 
{U — d Vi \l < i < n + 1} with this order is {U — d Vi |1 < i < n + 1}, and the 
initial ideal is generated by {U\l < i < n + 1}. Hence, the leading term of 
P — Q'kP'k *= P > ' w ^ n respect to the order -< must divide some U. However, 
the differential operator in D' which satisfies this condition is only 0. Then, we 
have P e D'F'. □ 

Corollary 1. The integration ideal J' is a holonomic ideal. 
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